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Observation of sub-Poisson photon statistics in the cavity-QED microlaser 
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We have measured the second-order correlation function of the cavity-QED microlaser output 
and observed a transition from photon bunching to antibunching with increasing average number 
of intracavity atoms. The observed correlation times and the transition from super- to sub-Poisson 
photon statistics can be well described by gain-loss feedback or enhanced/reduced restoring action 
against fluctuations in photon number in the context of a quantum microlaser theory and a pho- 
ton rate equation picture. However, the theory predicts a degree of antibunching several times 
larger than that observed, which may indicate the inadequacy of its treatment of atomic velocity 
distributions. 

PACS numbers: 42.50.-p, 42.55.-f 



Nonclassical light has attracted much attention in the 
context of overcoming the shot noise limit in preci- 
sion measurements and creating single photon pulses for 
quantum information processing 0. In quantum optics, 
one well-known source of antibunched light is in single- 
atom resonance fluorescence 0, 0, where antibunching 
occurs due to a "dead time" delay between photon emis- 
sion and atom re-excitation. The single-trapped-atom 
laser and similar setups for delivering photons on de- 
mand [ja, E[ exhibit photon antibunching essentially due 
to a similar process. The microlaser, on the other hand, 
generates nonclassical light via a very different process 
involving active stabilization of photon number, and re- 
markably, as shown below, photon antibunching and sub- 
Poisson statistics can occur even when the number of 
intracavity atoms greatly exceeds unity. 

The cavity-QED microlaser is a novel laser in which 
a interaction between the gain medium and optical cav- 
ity is coherent. Well-defined atom-cavity coupling and 
interaction time lead to unusual behavior such as multi- 
ple thresholds and bistability [g,|9|]- The microlaser has 
been predicted to be a source of nonclassical radiation 
due to an active stabilization of photon number at stable 
points. However, such predictions have generally been 
made on the basis of single-atom theory 12]. In this 
Letter we report the measurement of sub-Poisson pho- 
ton statistics in the microlaser even with the number of 
intracavity atoms as large as 500. 

The microlaser is the optical analogue of the micro- 
maser [ljj, in which sub-Poisson photon statistics has 
been inferred from the measurement of atom state statis- 
tics []JJ. The microlaser has the advantage of allowing 
direct measurement of statistical properties of its emitted 



field. 

Our experimental setup (Fig. ^) is similar to that of 
Refs. 0,111; 0- The optical resonator is a symmetric near- 
planar Fabry-Perot cavity (radius of curvature ro=10 cm, 
mirror separation L «0.94 mm, finesse F « 0.94 x 10 6 at 
A = 791 nm, cavity linewidth r c /27r=150 kHz). Barium 
atoms in a supersonic beam traverse the TEMqo cav- 
ity mode (mode waist u m = 41/xm) which is near res- 
onance with the ^o^Pi transition of 138 Ba (A=791.1 
nm, linewidth r a /27r «50 kHz). Shortly before entering 
the cavity mode, atoms pass through a focused pump 
beam which excites them to the 3 Pi state via an adia- 
batic inversion process similar to that described in Refs. 

In order to ensure coherent atom-cavity interaction the 
variations in atom-cavity coupling constant and inter- 
action time (or atomic velocity) have to be minimized. 
The sinusoidal spatial variation of the atom-cavity cou- 
pling constant g(r) along the cavity axis due to the cavity 
standing wave is eliminated by employing a tilted atomic 
beam configuration [lj, [l5j]. The remaining transverse 
variation of g is minimized by restricting atoms to the 
center of the Gaussian cavity mode (i.e., close to the 
plane containing the atomic beam direction and cavity 
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FIG. 1: Schematic of experimental setup. M: mirror, C: cav- 
ity mode, P: pump beam, 9: tilt angle, B: Ba atomic beam, 
Dl, D2: start and stop detectors, Counterl,2: counter/timing 
boards. 



axis) via a 250 /xmx25 fim rectangular aperture oriented 
parallel to the cavity axis. The resulting variation in 
peak coupling g is about 10%. The aperture is placed 3 
mm upstream of the cavity mode. The coupling at the 
center of the mode is 2go = 2n x 380 kHz. 

A supersonic beam oven similar to that of was 
employed to generate a narrow- velocity atomic beam. At 
highest temperatures it can produce a beam with a veloc- 
ity distribution of width Av/vq = 12% with v the most 
probable atom velocity and Av the width of the distribu- 
tion (FWHM). However, under these conditions the oven 
lifetimes were impracticably short. Instead, we used a 
lower temperature oven, which results in longer oven life- 
time but a broader velocity distribution, Av/v ~ 45% 
with vo ~ 750 m/s. The interaction time for an atom 
velocity at velocity i>o is t[ nt (= ^u m /vo) m 0.10 //s. 

In the tilted atomic beam configuration the microlaser 
exhibits two cavity resonances at u a ± kvo6, correspond- 
ing to two Doppler-shifted traveling- wave modes fl4lll5l]. 
The tilt angle 6 « 15 mrad corresponds to a separation of 
the two resonances by 2kvo9 ~30 MHz, and thus the con- 
dition for traveling- wave-like interaction pJI, kv^O ^> g, 
is satisfied. The cavity spacing is adjusted by a cylindri- 
cal piezo actuator to lock to one of the two resonances 
with a use of a locking laser. In the experiment, cav- 
ity locking alternates with microlaser operation and data 
collection. During data collection, the microlaser output 
passes through a beamsplitter and photons are detected 
by two avalanche photodiodes. 

The dynamics of the microlaser result from an oscilla- 
tory gain function associated with coherent atom-cavity 
interaction. For a two-level atom prepared in its ex- 
cited state and injected into a cavity, the ground state 
probability after the atom-cavity interaction time £i n t is 
given by P n sin 2 (\/n + lgtmt) , where P n is the intra- 
cavity photon number distribution function. When the 
mean number of photons (n) in the cavity is much larger 
than unity (i.e., semiclassical limit), as in the present 
study, the time variation of the mean photon number 
can be obtained by means of a semiclassical rate equa- 
tion E1II3 given by d(n )/dt = G((n}) - L((n}), where 
G(n) — (N) sin 2 [y/n + lgt int ) Aint the gain or emission 
rate of photons into the cavity mode and L(n) = T c n the 
loss with (N) the mean number of atoms in the cavity. 
The microlaser gain and loss are depicted in Fig. H{a). 
For comparison, G and L for a conventional laser are 
shown in Fig. H^b). 

Photon number stabilization or suppression of photon 
number fluctuations occurs when the gain has negative 
slope. Consider a momentary deviation in the cavity pho- 
ton number from a steady state value no- The gain and 
loss provide feedback, acting to compensate for the defi- 
ciency or excess of photons, and restore the photon num- 
ber to its steady state value in a characteristic time r c , 
the correlation time. The tendency to stabilize the pho- 
ton number is enhanced by the difference between the 
gain and loss as seen in Fig. [21 

Note that the rate to remove excessive photons or to 
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FIG. 2: Photon number stabilization. Solid line: gain G(n); 
dashed line: loss L(n). The restoring rate of the cavity-QED 
microlaser, shown in (a), for a momentary deviation An from 
a steady-state value no is d(L — G)/dn\ nQ . (a) Cavity QED 
microlaser has oscillatory G(n). (b) Conventional laser: G(n) 
approaches a constant value for large photon number and the 
restoring rate is dL/dn = T c . 



supplement deficient photons is not just T c as in the con- 
ventional laser, where the gain saturates to a constant 
value, but T c — dG/dn\ nQ > T c . The correlation time 

is then identified as r c = [d(L — G)/dn\ ~\ 1 . This en- 
hanced restoring rate is the source of suppression of pho- 
ton number fluctuations below the shot noise level and 
thus of the sub-Poisson photon statistics. In the semi- 
classical limit ((n) ^> 1), one can show that the Mandel 
Q parameter, defined as Q = (An) 2 / (n) — l with (An) 2 = 
(n 2 ) — (n) 2 the photon number variance |l8| . is approxi- 
mately given by Q ~ G f (jinU [T c — G'(no)] from the one- 
atom micromaser theory yj] with G'(no) = dG/dn\ nQ . 
Using the expression for r c above, we then obtain a sim- 
ple relation between the Mandel Q and the correlation 
time: Q = T c r c — 1. 

In the experiment we first measured (n) as the mean 
number of atoms (N) in the cavity was varied (Fig. Eta)). 
The photon number increases with (N) until it stabilizes 
(or saturates) around (N) « 200. Further increase in (N) 
results in a jump in (n). Similar jumps have been ob- 
served in micromaser experiments from sudden changes 
in atomic state 0] . The first direct observation of these 
jumps (or multiple thresholds) in the microlaser has re- 
cently been achieved @, 13 • 

The (n)- versus- (N) data can be well described by a 
quantum microlaser theory, i.e., the one-atom micro- 
maser theory |l2j extrapolated to large (TV), by the rea- 
sons to be discussed below. Fig. 0] shows the fit obtained 
by using this extrapolated quantum microlaser theory in- 
corporating atomic velocity distribution via averaging of 
G(n) over the velocity distribution. 

For r=0, the second order correlation function is re- 
lated to the photon number distribution for a stationary 
single mode by g^ 2 \0) = 1 + Q/(n). For our experi- 
mental parameters, sub-Poisson statistics requires sev- 
eral hundred photons to be present in the cavity. Thus, 
#( 2 )(0) is very close to 1, as Q > —1, requiring very 
low noise measurements of g^ 2 \r) for sub-Poisson statis- 
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FIG. 3: (a) Observed (n)-versus-(A^) curve. The solid curve is 
a fit based on the quantum microlaser theory, (b) Measured 
second-order correlation function g^ 2 \r) for (N)=12, 26, 68 
and 158. Each result is well fit by an exponentially decaying 
function. 



10- 



Experiment 
• Rate equation analysis 
- Quantum microlaser theory 



10- 




- Experiment 

Quantum microlaser theory 

T t from quantum theory 



50 100 150 200 250 

(N) 

(a) 



50 



n 1 r 

100 150 200 250 

(N) 



(b) 



FIG. 4: Microlaser correlation times (a) and Q values (b) 
versus (N), compared with theory. In (a), cavity decay time 
is represented by a horizontal dotted line. Solid line, quantum 
theory. Dashed line, semiclassical theory. In (b), solid line is 
the quantum theory. 



tics to be observed. To accomplish this, we developed 
a novel high-throughput multi- start multi-stop photon 
correlation system based on PC timing boards .20] and 
performed extensive averaging. With approximately 3 
MHz count rates on the two detectors and 300 sec to- 
tal acquisition time, the rms shot noise in g^ 2 \r) was 
0.00013. 

We measured g^ 2 \r) for seven representative points in 
the (n)-versus-(TV) curve. The results for four points la- 
beled A, B, C and D are shown in Fig. E^b). They are 
well fit by a function g^ 2 \r) = 1 + Coe -T / Tc , where neg- 
ative (positive) Co corresponds to antibunching (bunch- 
ing). From these fits we obtain the values of r c and Q 
shown in Fig. 01 Plots A and B in Fig. Ub), obtained 
in the initial threshold region, exhibit photon bunching. 
Data at C and D, from the region where photon number 
stabilization occurs, exhibit antibunching. The greatest 
degree of antibunching occurs at D, where (N) « 158 and 
Q = —0.13, corresponding to reduction in photon num- 
ber variance by 13% relative to a Poisson distribution. 

In Fig. 0fa), the observed correlation times are com- 
pared with the predictions by quantum and semiclassical 
theories. In the quantum theory, the correlation time is 
obtained via the quantum regression theorem [2l|. The 
predictions of the two theories are similar, and in agree- 
ment with experiment. For small atom number, r c is 
much larger than cavity decay time T" 1 but rapidly de- 
creases with increasing (N) to about half the cavity de- 
cay time. For the highest densities the r c may show a 
gradual increase. Although somewhat better agreement 
is obtained for the quantum theory, the observed corre- 
lation times are consistent with both theories, suggesting 
that the correlation time is primarily dependent on the 
dynamics of the mean photon number in the semiclassical 
limit. 

Fig. Elb) shows Q values for the different values of 



{N), along with the predictions of the quantum micro- 
laser theory, in which atomic velocity spread is included 
by integrating the gain function over the atomic velocity 
distribution function. The predictions of the theory do 
not agree with the data: the transition from super- to 
sub-Poisson distributions occurs at smaller (N) than the 
measured values, and the magnitudes of Q in the sub- 
Poisson region are about 5 times larger than those in the 
experimental results. By contrast, the predictions for 
(n) and r c versus (N) (Figs. Eta) and^fa), respectively) 
agree with the data. 

There may be several factors contributing to the dis- 
agreement in Fig.0{b). First, one may consider the cav- 
ity decay during the interaction time and simultaneous 
presence of many atoms in the cavity, both of which are 
not included in the quantum microlaser theory. However, 
we found from quantum trajectory simulations (QTS) 
|22l E3] that the inclusion of the cavity decay would in- 
crease Q by at most 0.1 0, |24|. We obtained similar 
results for the many-atom effects |2fj . 

Second, atomic velocity distributions may be inade- 
quately treated in the quantum theory by averaging of 
the gain G(n). We have performed additional QTS, 
which can more realistically describe velocity spread, and 
found that for large velocity widths QTS generally pre- 
dicts a broader photon number distribution than does 
quantum microlaser theory [2(| • The results suggest that 
the primary reason for the disagreement between theory 
and experiment is that the velocity distribution may not 
be treated adequately by the quantum microlaser theory. 
We are exploring modifications to the quantum theory to 
treat velocity distribution more accurately. 

Note that enhanced (reduced) restoring action in Fig. 
121 corresponding to a sub- (super-) Poisson distribution, 
occurs when dG/dn < (> 0), and thus results in a cor- 
relation time shorter (longer) than the cavity decay time. 
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In the semiclassical limit, Mandel Q is related to the cor- 
relation time by the relation Q = T c r c — 1 as discussed 
above. This relation is shown in Fig.|U where correlation 
times larger than the cavity decay time T" 1 correspond 
to Q > whereas the correlation time much shorter than 
the cavity decay time correspond to Q < 0. However, 
the observed correlation time when the transition from 
super- to sub-Poisson distributions occurs ((N) « 40) 
is shorter than the cavity decay time. This discrepancy 
might be due to the non-negligible atom/cavity damp- 
ing and significant atomic velocity spread, which could 
introduce additional fluctuations in the cavity field and 
thus a slightly enhanced restoring rate than the cavity 
decay rate would be needed in order to achieve a Poisson 
distribution for the cavity field. 

It may seem surprising that a single-atom theory can 
describe even the microlaser average photon number with 
a large number of atoms in the cavity mode. We have 
found that the cavity-QED microlaser can be well de- 
scribed by the (modified) single-atom micromaser theory 
p^ |. as long as gt int <C \J (n) |25| . which is well satis- 
fied in the present experiments. Under this condition, 
photon emission/ absorption by other atoms in the cav- 
ity does not affect the Rabi oscillation angle of a par- 
ticular atom interacting with the common cavity field 
since the angle change |A0| due to single photon emis- 



sion/absorption satisfies |A0| ~ gt- m t |An| /\/n + 1 <C 1 
for An = ±1. Therefore, the mean number of atoms (N) 
in the cavity becomes a pumping parameter in the frame- 
work of an extrapolated single-atom micromaser theory 

In conclusion, we have performed the first direct mea- 
surement of nonclassical photon statistics in the cavity- 
QED microlaser. The transition from super- to sub- 
Poisson photon statistics was observed as the mean num- 
ber of atoms in the cavity was increased. A minimum Q 
of —0.13 was observed for mean photon number about 
500. The observed correlation times and connection with 
the observed Q are consistent with a gain-loss feedback 
model. Values of Q reflect a lower reduction in pho- 
ton number variance compared to the predictions of the 
quantum theory; this disagreement will require further 
study. Our analysis suggests that in future experiments 
with a velocity distribution width of 15% it will be possi- 
ble to observe values of Q as low as —0.5. Other future di- 
rections include investigation of the microlaser field dur- 
ing jumps and measurement of microlaser lineshape |27| . 
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